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Abstract. Recently, the authors have proposed a new approach to the theory of 
random metrics, making an explicit link between probability measures on the space 
of metrics on a Kahler manifold and random matrix models. We consider simple 
pi I examples of such models and study the one and two-point functions of the metric. 

^ ■ These geometric correlation functions correspond to new interesting types of matrix 

model correlators. We provide in particular a detailed study of the Wishart model, 
where we determine the correlation functions explicitly. We find that the random 
measure in this model turns out to be concentrated on the background metric in the 



! large TV limit. 

oo" 
m 

^ . 1. Introduction 

This paper is a continuation of a series of papers by the authors on random Kahler 
metrics. In Refs. [9, 10] we have proposed a new method to define a random metric on 
^ ■ a Kahler manifold M. The purpose of this paper is to give a detailed study of simple 

[ models of random metrics where our methods apply. Our approach is based on the 

approximation of the Kahler metrics in a fixed Kahler class by the Bergman metrics. 
The latter are parameterized by positive-hermitian matrices Pij as follows 

9ab{z) = ^dadi,logSi{z)PijSj{z) . (l) 

Here the matrix P^j is contracted with a vector Si{z) (and its complex-conjugate Si{z)) 
defining a basis of holomorphic polynomials of degree k on the manifold, or more 
precisely, a basis of global holomorphic sections of some line bundle L'' on M. The 
dimension Nk of the vector depends on k and on the complex dimension n of the 
manifold. At large k, we have ^ k^ + 0{k'^~^). The positivity of the metric (1) is 
ensured by the condition that the matrix P is hermitian positive definite. Note that any 
such P can be written as P = A^A for some A G GL{Nk, C) defined up to multiplication 
by a unitary matrix on the left. Note also that constant rescalings of P do not change 
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the metric (1). Thus the space of Bergman metrics Bk is equivalent to the symmetric 
space 

Bk = SL{Nk,C)/SU{Nk). (2) 

The key fact [24, 25, 3] is that in the /c — > oo hmit, the space Bk of Bergman 
metrics of degree k tends in a strong metric sense to the infinite dimensional space of 
all Kahler metrics. As a consequence, many geometric properties of the Kahler metrics 
can be studied using the finite-dimensional spaces Bk, see e.g. [21] for review. With 
this correspondence in mind, we define a random metric on the Kahler manifold M as 
the large k limit of a random Bergman metric. Equivalently, formula (1) shows that 
defining a random Bergman metric is the same as defining a random positive definite 
hermitian matrix P. 

Probably the best known models of random metrics in physics are Polyakov's models 
of two dimensional CFTs coupled to gravity based on the Liouville action [22], see [8] for 
recent rigorous treatment. Other geometric actions, like the Mabuchi functional, were 
recently shown to play a role when non-conformal matter is coupled to gravity [9, 11]. 

Our goal in this paper is rather to study a much easier class of random measures on 
the space of metrics, which are natural from the point of view of the parametrization (1) 
of the metrics in terms of the matrices P. Namely, we consider the simple class of matrix 
models for which the random matrix measures diJ,{P) depend only on the eigenvalues of 
P. In other words, the measures we consider are invariant under conjugation of P by a 
unitary matrix U G U{Nk), 

dfi{U^PU) = di2{P) . (3) 

A particular example that we shall study in full details is the so-called Wishart model 
for which 

rf/i(P) = e^^*'^^[rfP]Haar, (4) 

where [(iP] Haar is the standard Haar measure on positive hermitian matrices defined in 
Section 3.1 and (7 > is a parameter. 

Matrix models of the type (3), and in particular (4), are standard matrix models 
whose large limits have been extensively studied using well-known matrix model 
technology. However, the geometric interpretation in terms of random metrics given 
by the formula (1) suggests to study an entirely new class of correlators involving the 
products of the logarithms 

J log Si{Zi)PijSj{Zi) ■ ■ ■ log Si{z^)PijSj{z„^) dfi{P) (5) 

depending on the points Zi,...,Zm on the manifold. These correlators are nontrivial 
because they do not depend only on the eigenvalues of the matrix P. If we decompose 
P = f/^Af/, where f/ is a unitary matrix and A is diagonal, we thus have to perform a 
nontrivial integral over U. This can be done using the Harish-Chandra-Itzykson-Zuber 
formula [13, 14]. The result depends on the coordinates Zp on the manifold via various 
inner products between the vectors of sections Si{zp). These inner products turn out to 
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be proportional to the Bergman kernel. We are then able to find the large k (equivalently, 
large A^^) limit of the correlators using the well-studied asymptotic expansions for this 
kernel [18, 16]. 

Our main results include explicit computations of one and two-point functions (5) 
for general measures of the type (3) for finite values of k, see Eqs. (41, 70) in the main 
text. The one-point functions turn out to be very simple and universal, independently of 
the precise form of the measure, see Eq. (44). The finite k two-point functions depend 
on the details of the measure, but their large k limit is universal and depend only on a 
model-dependent constant, see Eq. (76). In the case of the Wishart model (4), we go 
further and find an explicit formula Eq. (94) for the two-point function at finite k. The 
main feature of the measures considered here is that in the large k limit they turn out 
to be heavily concentrated around the background metric. In particular, the results for 
the Wishart model are consistent with the delta-function like distribution of the space of 
all Kahler metrics. This feature is an artefact of using the matrix P as a basic random 
variable. As was pointed out in [10], we expect much smoother limiting measures when 
using the "tangent space" random variable ^ log P instead of P. Nevertheless, we believe 
that the present work provides a necessary basis for further investigations of the random 
measures on Bk- 

The paper is organized as follows. In Section 2 we review basic facts on Bergman 
metrics and collect some formulas on the expansion of the Bergman kernel, which we 
use later in the text. We introduce the matrix model measures in Section 3.1. In Section 
3.2 we define the correlation functions and show how to integrate over the angular part 
of the measure. In Section 3.3 we compute the one-point function, and in Section 3.4 
we perform the angular integration in the two-point functions in the general case. For 
measures of the form 

rf/i(P) = e-*^^(^)[rfP]Haar, (6) 

where V is a. polynomial, we present the two-point function in a more explicit form 
in Section 3.5 using the method of orthogonal polynomials. In Section 4.1 we derive 
another representation for the two-point function, which allows us to derive the general 
form of its large k limit in Section 4.2. In Section 5 we focus on the Wishart model and 
derive the two-point function in an explicit form at finite k. We also compute its large 
k limit explicitly, checking our previous general results. 

2. Background 

2.1. Bergman metrics and Bergman kernel 

Here we collect some basic formulas from the setup of [10], which we will use in this 
paper. Consider a compact Kahler manifold M of complex dimension n. We fix a 
complex structure J on M and consider the set /C[cjo] of Kahler metrics in the cohomology 
class [uq] G H^'^{M, 27rZ), which can be parameterized by the space of Kahler potentials 

/C[,„] = {0 G C°°(M)/R ■.u^ = uo + idd(P > 0} (7) 
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Under the integrality assumption on [uq], there exists a holomorphic hne bundle L ^ M 
with Ci(L) = [ojq]. The space IC[^^] is thus identified with the space H of Hermitian 
metrics h = e~'^ho on L with the curvature = —ddlogh G [uq]. Consider then the 
kth power of the line bundle and the corresponding A^^- dimensional space H^{M, L^) 
of holomorphic sections. We choose a basis of sections {si{z)} = {si{z), sn^^z)), 
orthonormal with respect to the reference (background) metric on L^, 

^ s,{z)sj{z)h^o^^ = Sij, (8) 



V 



M 



where cuq is the Kahler metric on M corresponding to metric Hq the line bundle 

^oab = ^dadi\ogh'^. (9) 

The basis of sections {si{z)} is a vector in the linear space of dimension Nj. and it is 
defined up to a multiplication by a constant. Therefore it can be thought of as an 
embedding z — )■ [si{z), SJ^f^,{z)] of M into the projective space CP^*~^ of sections (the 
Kodaira embedding). All other choices of bases {si{z)} can be obtained by applying a 
linear transformation to the initial basis: Si{z) = AijSj{z), for A G GL{Nk,C). The set 
Bk of Bergman metrics is then defined as the puUback of the Fubini-Study metric^ from 

^<Ppab{z) = ^dadi\ogSi{z)PijSj{z), (10) 

where P = A'' A is a positive hermitian matrix in GL{Nk, C)/U{Nk). 

Note that a rescaling of P by a constant does not change the metric, so we can 
mod out by scalar matrices. Hence the space Bk of Bergman metrics is equivalent to 

Bk = SL{Nk,C)/SU{Nk). (11) 

The key fact, shown in [24, 25, 3], is that for large k the space of Bergman metrics 
approximates the space }C[uio] of all Kahler metrics 

/C[,„] = lim Bk. (12) 

k—¥oo 

The proof of (12) is based on the asymptotic expansion of the Bergman kernel on the 
diagonal [25, 3, 15, 7]. For the orthonormal basis of sections as in Eq. (8) the latter has 
the form 



Pk[Z) 



\s{z)\X ■■= ^*(^)^i(^)^o = ^"(1 + i^(^o) + 0(1/^)). (13) 



2k 
=1 



Here all the terms that appear on the right are various curvature invariants, depending 
only on the metric uq. The total number of sections Nk is given by the integral of the 
Bergman kernel over the manifold. It grows at large k as 

Nk = A;" + ^ci(M)r"i + 0(r-2) (14) 

see [15] for the subleading terms. 

I With some abuse of notation we will make no distinction between the metric g and the corresponding 
Kahler form uj, connected as io^i ~ Wat- 



Simple matrix models for random Bergman metrics 5 
2.2. Off-diagonal Bergman kernel 

Another important function, which will appear later is the normalized off-diagonal 
Bergman kernel 

\{s{zi),s{z2))\'^ 1 \^ - f \ / \- r \ r \ nK\ 

Note, that this is nothing but the squared norm of the inner product of the two vectors of 
sections Si{zi) and Si{z2)- If ^12 is the angle between the two vectors, then the Bergman 
kernel is given by 

2n ( ^ |(g(^l),g(^2))P 

cos ^12^1,^2 = 16 
\s{Zi)\^\s{Z2)\^ 

We will need the Taylor expansion of (16) for small distances between the points Zi, Z2, 
which can be obtained following [18, 16, 23]. First, we rewrite (16) in terms of the 
Calabi's diastatic function, defined as follows. Any Kahler metric can be written in 
local coordinates as 00^1 = dadi^{z, z), for some potential ^{z,z). Consider ^{z,z) 
formally as a function of two independent variables z and z. Then the Calabi's diastatic 
function is defined as 

D{zi, Z2) = ^zu zi) + Hz2, Z2) - $(zi, Z2) - $(22, zi). (17) 

This function is an invariant of the Kahler metric, independent of the choice of local 
coordinates [2]. Consider now the Bergman metric 

u^^ = ^dd\og\s{z)\\ (18) 

corresponding to the identity matrix P = I (10). The normalized Bergman kernel (16) 
can then be written as 

cos^ ^12(^1, Z2) = e-^^'^''^'^\ (19) 

where 

Di{z,, Z2) = ^(log |.(zi) r + log \s{z2)\' - log I {s{z,), s{z2m (20) 

is the diastatic function defined with respect to the metric (18). This function equals 
zero only for coincident points Zi, Z2- In normal coordinates around a point = we 
have the expansion 

$(z, z) = \z\' + Ji?„,,5(0)^'^^V# + 0{\z\') (21) 

With some abuse of notation we assume that zi and Z2 are the coordinates of the points 
in the chosen local chart. Therefore we get the following Taylor expansion for small 
distances between zi and Z2 

Diizu Z2) = \z,-z2\^+\r,-,,m {4444 + 444^2 - 4444 - zi442^^+- ■ ■ (22) 

where all the quantities here are defined with respect to the metric (18). The expansion 
of the off-diagonal Bergman kernel follows immediately from the above formula. 
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3. Eigenvalue measures on Bergman metrics 
3.1. General form of the measure 

Many aspects of the infinite-dimensional geometry of /C[^(,] (7) can be approximated 
by the symmetric space geometry of Bk (11). This is the main idea behind the Yau- 
Tian-Donaldson program in Kahler geometry [24, 6, 21]. In Ref. [10] we proposed to 
define probability measures d/i^ on random Kahler metrics as large k limit of probability 
measures on the finite-dimensional symmetric spaces Bk, 

/ F{u^)diJ,^ := \im / Fk{uj^p)dnB^{P). (23) 

Here F{u^) schematically denotes operators, whose correlation functions we want to 
compute. Some examples of physically interesting measures dn^ are described in [10]. 

In particular, the natural metric on lC[uio\i the Mabuchi-Semmes-Donaldson metric, 
arises in the large k limit from the invariant metric on Bk [1], see [10] for the details. 
The volume element of the invariant metric is the Haar measure 

WW. = [dP] (24) 

f^^^ " fdetPFfe ^^^^ ' ' ' ^^^^^^ n ^ ^ ' 

which we use as a basic building block for the random measure on the space of Bergman 
metrics. Recall for future reference the standard eigenvalue decomposition of the 
measure 

[dP] = A\\) [dX][dU], A(A) = J](A, - A,), [d\] = l[dK, (25) 

i<j i=l 

where P = WAU, [dU] is the Haar measure on U{Nk) and the eigenvalues are positive 
definite A^ > 0. 

Simplest random matrix measures on positive hermitian matrices are the eigenvalue 
type measures. These are invariant under the U (Nk) conjugation P — )■ U'^PU. Therefore 
they have the general form of the product of the Haar measure (24) and an eigenvalue- 
dependent weight 

rf/i(P)=^(A)[dP]Haar, (26) 

where we use the standard decomposition of the hermitian matrix P = U'^AU into 
eigenvalue and angular part. We will also assume that J^{X) is a completely symmetric 
function of the eigenvalues. The measure on Bk can be obtained by gauge-fixing the 
scale invariance of Bergman metrics under P ^ c - P. This can be done by constraining 
the measure (26) to matrices P with unit determinant. We thus in practice will consider 
the gauge-fixed measures dfiBi,{P), defined as 

dfiB,{P) :=^B,(A)[rfP]Haar, (27) 
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where the new eigenvalue weight J-jg^. on Bk equals J-' multiplied by the delta-function 
constraint 

J-s,(A):=^(A)5|^^logA.j. (28) 

Note, that resulting measures are still strictly speaking of the form (26). The 
normalization is always chosen according to 

/ rf/iB,(P) = l. (29) 

Another useful parametrization of the Bergman metrics is in terms of GL{Nk, C) matrix 

y4., clS 

P = A^A, (30) 

with residual symmetry A — > VA, where V G U{Nk). The measure [dP] can be 
represented as [dP] = [dA]/[dV], where [dA] is the Lebesgue measure on complex 
matrices 

[dA] = JJ Re Aij dlmAij, (31) 

l<i,j<Nk 

and [dV] is the Haar measure on U{Nk). 



3.2. Correlation functions 

The basic set of correlation functions consists of the products of Bergman metrics at 
different points 

Efc ^<j>paibr (^l) • • • ^(t>pambm (^m) : = 

= da^Bl^ . . . da^Bi^ I (f)p{Zi)...(j)p{Zm)dflB^{P) (32) 

where we introduced the function 

(f>p{z) = ^\ogSi{z)PijSj{z), (33) 

called the Bergman potential. It would also be interesting to consider more complicated 
geometric correlation functions, e.g. involving the Ricci tensor R^i = —dadi,logu^. 
Following [9, 10], the main idea is that in the large (equivalently, large k) limit 
matrix model correlation functions on Bk will tend to exact correlation functions for 
corresponding limiting random measures on the full space of Kahler metrics /C[a;o] 

Hai6i(2;i) • ..u^^^i^iz^)) := lim^EkU^^^^i^izi) . . .u^^^^i^{z^). (34) 

For the eigenvalue models correlation functions depend explicitly on a choice of 
basis of sections, and therefore on the background metric cuq (8), corresponding to that 
choice. Note, that dependence of the Bergman potential (33) on the vector Si breaks 
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the U{Nk) invariance of the integral, making the integration over the angular variables 
nontrivial. Luckily, for the models of the type (26) the integration over the angular part 
can be performed explicitly, using the Harish-Chandra-Itzykson-Zuber integral [13, 14]. 
Indeed, we can use the following integral representation of logarithm 



lo. 



1 C"^ 

,ga = - + 7-/ x*^ie~"^rfx + C(t), (35) 
^ Jo 

for small t, and 7 is the Euler constant. Since the first two terms here are independent 
of Zi, we can express the correlator (32) as follows 

Efc W^paiSi (^1) • • • ^c^pa^b^ i^m) = (36) 

\J^daA...da^d,, \im II dx.x^' e-Sr^^^^^"(^')^^(^^)rf/ie,(P). 



Assuming the ansatz (26) and using the eigenvalue decomposition (25), the matrix 
integral in Eq. (36) reads 



f e-^i=i"''^"("')^^("')rf/iB,(^) 

g-tr UAWir-=r-M-s{z,))jr^^^x)A\X) [dX][dU], (37) 



and the unitary integral here is of the Harish-Chandra-Itzykson-Zuber type 



U{N) Volf/(iV) J A(a)A(6) 

Here and bj are the eigenvalues of some hermitian matrices A and B. Therefore, 
the angular integral in Eq. (37) depends on Aj and on the eigenvalues of the following 
matrix 

m 

^(^\^ = J2xM^g)sj{z,). (39) 

q=l 

This matrix is degenerate and has in general m non-zero eigenvalues. 
3.3. One-point function 

The formula (38) simplifies further for small m, i.e. for small number of eigenvalues of 
(39). In particular, the expectation value of a single metric can be found exactly 
for all eigenvalue type measures (26). The matrix has a single non-zero eigenvalue, 
which is equal to x|s(2;)p, and the angular integral in (37) gives 

where Aj(A) = Y[j<i jij^i^'^j " -^O Vandermonde determinant with omitted z'th 

eigenvalue. We get for the expectation value of a single Bergman metric 



^koo^jz) = --ddYim / dxx*-^ / e-"^"(")^'*(")ci/XB,(P) = -T^^lim 



\s(z)\-'' 
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/ dxx'-"^^ I ^(-l)-%(A)e-^^'J-e,(A)A(A)[rfA].Volf/(iV,) 
Jo Jm^^ 

1 - 

= --dd\im\s(z)\~^\Nk - m-lf'^-^CN, x'-^^e-^'dx (41) 
k t^o 

= CN^^-ddlog \s{z)\^ = cn^uj4,j{z), 

where the Bergman metric u^^ corresponds to the identity matrix P = I (18). We also 
sometimes omit the space indices, we it cannot lead to confusion. The constant cn^ 
equals one due to the normalization condition (29) for the measure 

p Nk „ 

= / . 5^(-l)^'=+'A,(A)Af-^J-^,(A)A(A) [dX] ■ VolUiN,) = / rf/i^, = 1, (42) 

since the sum over eigenvalues here is just the expansion of the Vandermonde 
determinant along a column. Finally, using the asymptotic expansion (13), we get 

u^^iz) = uoiz) + ^ddhgpkiz) = uoiz) + Oil/k^), (43) 

where Uq is the background metric (8). Therefore, in the large k limit the one point 
function of the Bergman metric is equal to the background metric 

{oj^{z)) = lim EkU^p{z) = Uo{z). (44) 

k—>oo 

This result is universal for models of the type (26). However, the case of the two-point 
function is model dependent. In the remainder of this section we analyze the two-point 
function for generic eigenvalue models, and in the Section 5 we derive it explicitly in 
the case of Wishart measure. 

3.4- Two-point function 

For the two-point function the matrix (39) has two non-zero eigenvalues 

01,2 = ^(a;i|s(2i)p + X2|s(22)n± 

^^J{xMzlW-^2\s{z2Wy + M{s{zl),s{z2))fx,X2. (45) 

Performing the angular integral in (37) and rescaling xi^2 2:1^2 1 ■5(2:1,2) |~^, we get for 
the two-point function 

EkUJ^p{Zi)uJ^p{z2) = 

-| poo poo 

= -ddUMz, lira \s{z,)\-^'^\siz2)r^'' dx.x'r' dx^x'r' 
k^ *i.*2^o Jo Jq 

4 ^W^n(X)e-*-*'-T,JXm ■ yoW(N,), (46) 

where A12 is the Vandermonde determinant for eigenvalues A3, A4, . . . Aat^., and we 
replaced (45) by the rescaled eigenvalues 

^(Xi + X2 ± {Xi - X2y + 4 C0S2 612X1X2). (47) 



± 
2 
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Now we can see that all nontrivial coordinate dependence of the two-point function 
is contained in the norms and in the off-diagonal Bergman kernel kernel 

cos'^ 9i2{zi, Z2) (16). From the expression (46) we can read off the general structure 
of the two-point function. The small t expansion of the xi, X2-integrals in Eq. (46) 
starts in general with singular terms 

+ + + /fc(cos2 ^12(^1, Z2)) + 0{tiM). (48) 

tlt2 ti 12 

and the coefficients in front of the singular terms are independent of the coordinates on 
the manifold. Therefore the two-point function has schematically the following form 

EkUj^p{zi)u^p{z2) = uj^j{zi)uj^j{z2) + ^dd\,,dd\;,Jk{cos^ Ouizi, Z2)), (49) 
and the scaling behavior of the function fk is to be determined. 



3.5. Single-trace potentials 

The two-point function (46) is simplified further when the eigenvalue weight in Eq. (26) 
has the form of the exponent of a single-trace potential 

jr(A) =e-*'"^(^) (50) 

The expression (46) can be simplified using the method of orthogonal polynomials [19, 5]. 
Suppose {pn{\)}, n = 0, Nk — I is the basis of orthogonal polynomials satisfying 

p„(A)p„(A)e-^(")A"rfA = h^{a)Sr^m (51) 







and normalized according to p„(A) = A*^ + . . . The parameter a here is assumed to be 
pure imaginary and is used for the gauge-fixing. Indeed, we can connect the measures 
(26) and (27) as 

/ioo 
da{deiPYd^i{P). (52) 
-ioo 

Then, using the relation 

A(A) = detP„_i(A„), (53) 
the normalization condition (29) reads 

/ dfiB,{P) = Nkl / da II hn{a) ■ Yo\U{Nk) = 1. (54) 

JBk J -ice, 

Let us now return to Eq. (46) and rewrite the Vandermonde determinants in the 
integrand using Eq. (53). 

The orthogonality condition Eq. (51) then allows us to perform all but two 
eigenvalue integrals, and we get 

Ek^<t>p{Zl)i^^p{z2) = 
T roo poo 

-ddUddU, lim \s{z^)\-^'^\s{z2)r''' / dx^x'r' / dxs^"' 
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nV'i - ^2) i-ioo 
Jo Jo 

■ (p^,~2(Ai)p7Vfc-i(A2) - m-i(^i)m-2(A2)). (55) 

We will use this representation later in order to derive the two-point function in the 
case of the Wishart potential. 

4. Integrating out matrix subspaces 

4.1. One and two-point functions 

Here we use the fact, that the matrix in Eq. (39) has only m eigenvalues, in 
order to obtain another representation of the two-point function, useful for the large k 
analysis. The idea is to integrate out subspaces of Bk, corresponding to the directions 
where has zero eigenvalues. The resulting formula will allow us to compute the 
general form of the large k limit for measures of the type (26). 

In order to illustrate the idea, let us start first with the one-point function, and use 
explicitly the decomposition of the matrix P into eigenvalues and angular variables 



EkUJ^p{z) = jdd [ log Si{z)UlXjUjiSi{z)dfiB^{P). 
Jb. 



(56) 



Consider now the unit vector Si{z)/\s{z)\ G and let G U{Nk) be any unitary 
matrix such that {U'^eiji = Si{z)/\s{z)\ G C^*, where ei is the unit vector pointing in 
the direction i = 1. We get 

EkU^^iz) = u^^iz) + ydd [ \og{UU'ei)]Xj{UU''ei)jdfis,{P) 

JBk 

= u^j{z) + -dd / \og{Uei)\\j{Uei)j dnB,{P) = w<^,(^), 
^ Jb, 

where in the second line we use the [/(A'^fc)-invariance of the measure, and the resulting 
integral is independent of z. This is a short way to reproduce the result (41). 
Now we would like to calculate the two-point function 

p{zi)(l)p{z2) = J \ogSi{zi)UljXjUjiSi{zi)\ogSi{z2)UljXjUjiSi{z2)diJ,^^^ (57) 
using the same idea. We can find a matrix U^'^'^'^ such that 
[U ei i - — — -, 

(f/^i'^^(ei cos ^{zu Z2) + 62 sin ^{z^, Z2))), = (58) 

\S{Z2)\ 

and taking the inner product of these two vectors we immediately get 

COS ip[Zi,Z2) = I . .,^1 . = COS 6i2{Zi,Z2) (59) 

\s(zi)\^\s(z2)r 



E 



k ■ 
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the off-diagonal Bergman kernel (16). Note that there may exist many such matrices 
f/zi,22_ Using the U{Nk) invariance of the measure, the integral (57) simplifies 

1 r 

EkMzi)M^2) = 0/(^1)0/(^2) + / log VA,|(t/t/^^'^^e0.r 
•log^A,|(t/?7^^'^Hei cos^i2 + e2sinM),|2rf/iB,.(^) = 0/(^i)0/(^2) + 

i=l 

/ log VA,|(f/ei)i|2log VA,|(f/(eicosei2 + e2sin^i2)).|'rf/iB,(P). 

■^'3^ i=l i=l 

The last line can be expressed in terms of the matrix elements of P 
Efc 0p(zi)0p(2;2) = 0/(^1)0/(2^2) + (60) 
/ log Pn log(Pii cos^ 6^2 + {Pi2 + ^2i) COS 012 sin 9^ + P22 sin^ O^) dfie^ {P) 

The integrand here depends only on the 2x2 block 

P(2) := I I e GL(2,C)/f/(2) =: P2. (61) 

\ ^21 -'22 / 

One can check that P^^) is a positive Hermitian matrix since for any v G C^, 
(p(2)^;^ v) = (p(^,, 0), {v, 0)) where {v, 0) G C^'^ is a vector whose last Nk — 2 entries are 
zero. Thus, we have the map 

U2 : GL{Nk, C)/U{Nk) ^ GL{2, C)/f/(2) (62) 

which takes the upper left 2x2 block from P. Indeed, if 112 '■ C'^'-' — is the map 
n2(2;i, Z2, z') = {zi,Z2, 0), then P^^^ = n2Pn2 and it is clear that this is a positive semi- 
definite Hermitian operator on C^* whose restriction to the range of 112 is positive. We 
note also that g embeds in the closure jj^ of positive semi-definite Hermitian matrices 
if we put the zero matrix on the lower diagonal Nk — 2 block. 

Since the integrand in Eq. (61) is manifestly a function of the three real variables 
Pn, RePi2, P22, we have 

/ log Pn log(Pn cos^ 612 + (P12 + P21) cos Ou sin 612 + P22 sin^ Ou) rf/ig, (P) (63) 

= I log Pn log(Pn cos^ 612 + (P12 + P21) cos ^12 sin 6^2 + P22 sin^ 612) d^ir,,k{P^^^). 
where dfx-p^^kiP^'^^) is the pushforward measure on V2 under the projection (62), 

f f dfJ^Bk jP) \ rjp(2)n 



df^v^AP''^) = ) [dPn (64) 

where the integral on the right runs over all matrix elements of P except those of V2 
and [(iP^^^] is the standard volume element on V2, given by 

[rfp(2)] = dPiidP22dRePi2dlmPi2- (65) 
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The t/(2)-invariance of this volume element immediately follows from (25). We claim 
that dfip^^kiP^^'^) is also f/(2) invariant under P^^) ^ f;{2)tp(2)f/(2)_ jn^jeed, we may 
embed U(2) U{Nk) as the upper 2x2 block, with the identity matrix on 
the lower diagonal. The original measure dfi^^^P) is U{Nk) invariant and therefore is 
invariant under this subgroup. Then this fact and U{2) invariance of (65) imply U{2) 
invariance of d^p^^kiP^'^^)- 

We may write the integral (63) in the eigenvalue and polar coordinates as follows. 
We write P^^) = f/(2)t^(2)(^^)f/(2) (non-uniquely) , so that f/(2) e U{2),D^^\X) = 
diag(Ai, A2), ei = (1, 0), 62 = (0, 1). Then the measure (64) has the form 

d/ip,,fc(p(2)) = J^2,kWdXid\2[dU^% (66) 

and [(if/^^^] is the Haar measure on U(2). The function J^2,fc(A) here contains the 
information about the original measure on Bk- Now we can rewrite the integral (63) as 

logPn log(Pn cos^ 9^ + (P12 + P21) cos ^12 sin ^12 + P22 sin' ^12) rf/ip„fc(P('^) = 

log 1^(2) (^)f;(2)g^ |2 1^(2) ^^-^^(2) ^ ^-^ |2 [dU^^)]J^2,kW>^ld^2. 

Consider first the inner integral here 

/2,fc(A,cos2^i2) = (67) 
log \D^^\\)U^^'^e^\Hog\D^^\\)U^^\e^cosei2 + G2^i^0i2)?[dU^^\ 

51/(2) 

Note that scalar matrices in U (2) commute with diagonal matrices, and clearly cancel 
in the integrand, so the integral over f/(2) is reduced to the one over SU{2). We can 
parameterize the latter as 

U = a^'kpa^, 

where 

' e'"^ \ , _ / cos /3 sin (3 
e-^"^ J ' ^ ~ i - sin /3 cos /3 

and the Haar measure on SU (2) equals 

[rff/(2)] = A-d4>' d4> sin/3 0<(f),4>' <2n < /3 < tt. (68) 

Since commutes with D^'^\\), the integrand in (67) is independent of The first 
logarithm is also independent of a^, and we can simplify the second logarithm in order 
to obtain the following integral, now over 

/2,fc(A, cos 6*12) = 

= / log |P>(2)(A)A;^ei|2log|P>(2)(A)V0(eicos^i2 + 62 sin^i2)|'# sin/3 

= / log |P)(2)(A)A;«ei|2log|P)(2)(A)A;«(eie^'^cosei2 + e2e-''^sin^i2)|'rf(/) sin/3d/3. 

4:71 Js2 
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Using the Jensen formula 



a + _ 52 

log(a + fecos</))(i</) = vrlog , a > |6| > 



we can perform the 0-integration 

/ log \D^^\\)kp{eie'^ cos + esc"^^ sin ^12) = 
Jo 

A + _ 52 



/ log(y4 + 5cos2(/))d(/) = 27rlog. 



2 

where we introduced the following notations 

A= \D''^\X)kpei\\osHi2 + \D'-^\X)ki3e2\^sm^ei2 = 
= (A? cos^ P + Xl sin^ (3) cos^ ^12 + (A? sin^ f3 + Xl cos^ /3) sin^ 6^2, 
B = 2{D^^\X)k(sei,D^'^\X)kpe2)cos9i2 sin^i2 = 2(A? - A^)cos^i2 sin^i2 cos/3sin/3. 
This removes one integral and leaves us with 

1 /-T A \ J A2 __ D2 

J2,fc(A, cos^ ^12) = - / log(A? cos^ /3 + A^ sin^ /3) log sin /3 (69) 

2 Jo ^ 

Our final answer for the two-point function (57) is given by 

Efc 0p(2l)0p(2;2) = h{zMl{^2) + i / /2,fc(A, COS^ ^12)J^2,fc(A)c/AirfA2, (70) 

<^ Jul 

where the function J^2,k is defined in (66). Now we are ready to perform the analysis of 
the large k of the two-point function. 



4-2. Large k limit in general case 

Representation (70) allows us to study the large k limit of the two-point correlation 
function for measures of the type (26). Note that the main coordinate dependence of 
the two-point function is encoded in the normalized off-diagonal Bergman kernel, which 
depends exponentially on k (19). Therefore we distinguish the following two cases when 
considering the large k limit: 

(1) kDj{zi,Z2) — !■ 00, or equivalently cos^6'i2 — )■ as A; — cxd, 

(2) kDi{zi, Z2) — !■ 0, or equivalently cos^ 6^12 — )■ 1 as A; — )■ 00. 

In the case (1) the distance between the points Zi and Z2 is finite and does not scale 
with k, so we can use C, = cos^ 612 as a small expansion parameter. In this case we obtain 
a limiting expression outside the diagonal Zi = Z2- In the case (2) the distance between 
zi and Z2 decreases effectively as l/A;^/^ or faster. Therefore, in this case e = kDj{zi, Z2) 
becomes the small parameter. This limit will produce contact terms, supported on the 
diagonal ^1 = ^2- 

The key fact is that the integral (69) is well-behaved and admits a Taylor expansion 
both around = cos^ ^12 = and e = kDj{zi, Z2) = 0. Our only assumption here is 
that this Taylor expansion is allowed to be used in the integral over eigenvalues in Eq. 
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(70). Therefore, in order to study the large k hmit in the case (1) we can expand the 
second logarithm in (69) around C = 

l2,k{X,cos'eu) = /2,;c(A,0) + A/2,(A,C)|^^oe-^^^(^-^^) + ^(e-^^^^^^-^^)), (71) 

Where /2,a;(A, 0) and the first derivative at zero I2 i.{X, 0) are given by certain convergent 
integrals, whose explicit form we do not need for this analysis. For finite Dj{zi,Z2) 
all the terms starting from the second on the right go to zero uniformly as /c — )■ 00. 
Therefore, we conclude that 

{uj^{zi)uj^{z2)) = lim EkUp{zi)up{z2) = ^0(^1)^0(^2), Di{zi,Z2) > 0. (72) 

fc— >oo 

In the case (2) we use kDj{zi, Z2) as the small parameter. Since Di ~ \zi — 22P + 
0{\zi — Z2\'^)-, as follows from Eq. (22), and since the correlation function of two Bergman 
metrics contains four 2;-derivatives, it suffices to expand (69) up to {kDjY order 

= ^cj>iarhi^l)^<t>ia2hi^2) + ^ lim da,di,^da2\ / ( -^l2,k{\, e-')\^^Q kD j{zi, Z2) + 

+\^,h,k{\ e~') l^^o {kDi{zi, Z2)f + 0{{kDi{zi, Z2)f)] T2,k{\)d\id\2 = 



2de^ '''' ' yi,=o ^ 

= ^(/./ a 1 51(2^1) 0262 (^2) + Ck{u;^jaibi^<l>ia2b2 + ^<t>ia2bi^(t>iaib2)\zi^zu^2 (73) 

where 5^1,22 equals one on the diagonal and zero outside. The constant Ck is given by 

r d^ 

Cfc= / T^^2,fc(A,e-^)|^^oJ-2,fc(A)rfAidA2, (74) 

and the final result depends on its asymptotic value 

c = lim Cfc. (75) 

k—^oo 

Combining the on-diagonal (73) and off-diagonal (72) results and taking the large k 
limit we get 

{<^(l>aiai{^l)<^(f)a2a2{^2)) = Efc OJ^pQ^Q^ (^i )0L;(^pa2a2 (^2) = 

k—^00 

= i^Oaibi (^1)1^00252 (^2) + c{uQaibi^0a2b2 + '^0a2bi '^00162 ) Lj '^^i ,Z2 • 

(76) 

which is the general form of the two-point correlation function for the eigenvalue type 
measures (26). 

Let us also note that in addition to the off and on-diagonal limits (1) and (2), one 
can also consider the scaling limits of the correlation function (70). They correspond 
to the "near-diagonal" asymptotics of Ref. [23] and reveal a fine structure behind the 
contact terms in (76). The idea is to make the distance scale as D'j{zi,Z2 ) < bHogk/k 
for some constant 6 as A; — oo. In local coordinates around zq this amount to taking 
Zi = Zq + u/y/k and Z2 = Zq + v/ y/k, where |f | < by/\og k. The Bergman kernel (19) 
then behaves as 

cos' ^12(^0 + u/Vk, zo + v/Vk) ~ e"!"-''!', (77) 

and the Taylor expansions of the previous cases no longer apply. Existence of such 
scaling limits is an intrinsic feature of the models of eigenvalue type. 
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5. Wishart ensemble of random Bergman metrics 

5.1. Wishart distribution 

Now we would like to consider an explicit solvable example where correlation functions 
can be computed explicitly for finite k. 

Consider the Wishart distribution on positive hermitian x A^^ matrices, given 

by 

d/i„(P) = ^e'^'^^'idetPridP], (78) 

and depending on the two parameters a and g. It can be equivalently written as in 
Eq. (30) in terms of the GL[Nk., C) matrix A, taking P = A^A with residual symmetry 
A^VA 

dfiM^A) = ^e-^^^^'^{detA^Ar^^, (79) 

Zg [dV\ 

where the measure [dA] is defined in Eq. (31). The constrained Wishart measure (27) 
on Bk can be obtained from (78) by the Laplace transform 

/ioo -| 
dfXa{P) da = — *''^5(logdet P)[dP] (80) 

and we are left with a single parameter g. The constant Zg is determined by the 
normalization condition (29), which reads 



Zg= / e~^''^6{logdetP)[dP] (81) 
More explicitly 

/ioo f N 

da g-^^(^^+'^) / e- I A(A) |2 n A^dA, ■ Vol U{Nk) = 
-ioo Jr^" fj[ 

/ioo Nk 
rf^^-A^.{iv.+a) -Qr^^. ^ ^ ^) . voif/(iv,) = 
■ioo 



-N? 

g 



k 



i=l 



<-°(^7^''|l,2,...,iV,).Volt/(iV,), (82) 



where the Meijer G-function in the last line is defined by the integral in the second line. 
The integral in the second line is carried out using the Selberg formula, see e.g. Eq. 
(17.6.5) in Ref. [19]. We also define another normalization constant Zga as 

Nk 

Zg,a = n + l)r(j + «) ■ Vol U{m) (83) 



Simple matrix models for random Bergman metrics 17 
5.2. One-point function 

Correlation functions in the Wishart model can be computed using the GL{Nk,C) 
representation of the measure (79). Let us consider first the case of the one-point 
function. We introduce the notation for the expectation value with respect to the 
unconstrained measure (78), related to the expectation value with respect to the 
constrained measure (80) via the Laplace transform 



Efc...= / daEl... (84) 

J —ioo 

Then for the one-point function we get 

1 - /"°° 

Eloo^^ = --dd\im dxx'--'Ele~'''^'^^'^'\ (85) 

Since the expression in the exponent is linear in P we can rewrite the expectation value 
on the right as 

^a^-xs{z)Ps{z) = ^ f e-^*"^^^'+f"^"(^))(detP)'^ [dP] = 

f^-9trAUil+^s.H^))^^^^^tAyM, (86) 

Zg J [dV\ 

where in the second line we use the GL{Nk, C) representation (79). The positive definite 
hermitian matrix in the exponent of Eq. (86) can be written as 

lij + ^Si ■ Sj = {BB% (87) 

for some matrix B G GL{Nk, C). 

Using the multiplication covariance property [^(SA)] = {det B^ B)^'' [dA], we 
immediately get 



E«w. = - ^^aiim / dxx'~^[det(I +-s-s)y 
k Zg t^oJo g 



iNk+a) 



i^^aiim r dxx'-Ul + -\s\^r^'''^''^dx 



k Zg t^oji^ g 
= -l^dd\im\s\~^'g'B(Nk + a,t) = r^dd\og\s\^ (88) 

k Zg k Zg 

where B[Nk + a,t) is the beta-function, with the following behavior near t = 

, T(Nk + a-t)T(t) 1 , 

The correlator (84) on Bk agrees with the previous result (41), 

/ioo 2 
daElu^^iz) = -dd\og\s\^ = cuoiz) + Oil/k^). (89) 
-ioo 
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5. 3. Two-point function 

For the two-point function the calculation is analagous. We write 

Etu^^{z,)u^^{z2)= (90) 

-| POO POO 

The matrix BB^ (87) in this case equals 

BB^ = I + ^s{z,) ■ 7s{z,) + ^s{z2) ■ rs{z2) (91) 
9 9 

and the expectation value of the exponent reads 

■^a ^-^siPs{zi)-^S2Ps{z2) _ 

= %^[det(/ + ^s{z,) ■ s{z^) + ^s{z,) ■ si^^W^""'^^^ = 

Zg 9 9 

= ^(1 + -Hzi)\' + ^Hz2)\' + ^{Hz,)\M^2)\' - |(s(.i), siz.mr^'^^^^K 

Zg 9 9 9^ 

Rescaling xi_2 a^i,2|s(2i^2)|~^/fi' we reduce the two-point function to 

EfcW</,p(2;i)w0p(z2) = 

- l^OdUddU \im\s{z^)\-''^g'^\s{z,)\-''-g'^ 



-oo 

-{Nk+a) 



poo poo 

■ / dxix{^^^ / (ix2 X2^~^(l + Xi + X2 + X1Z2 sin^ ^12) 
Jo Jo 



(92) 



We also derive this formula in the Appendix in a more direct way, using the method of 
orthogonal polynomials. 

Now the integrals over xi,X2 in (92) can be performed 

Elu^^{z,)u^^{z2)= (93) 
- ^ ^''""ddUddU. lim \siz,)\-''^g'^\s{z2)\-'''g''- 



B(ti, Nk + a- ti)B(t2, Nk + a- t^) 2^1(^1, ^2, Nk + a; cos^ ^12)- 

In order to compute the small t limit in Eq. (93) we use the following representation of 
the hypergeometric function 

B(t2, iVfc + a - ta) 2Fi{ti,t2, Nk + a; cos^ 612) = 

= I dxx*2-i(l - x)^'=+''-*^~^(l -xcos2ei2)"*S 
Jo 

see e.g. Ref. [1]. Then for the expression in Eq. (93) we get 
lim \s{zi)\-^''g''\s{z2)\-''''g'' B(ti, N^ + a- ti)B(t2, N^ + a- ^2) 

tl,t2-!>0 

■2Fi(ti,t2,Nk + a; cos^6'i2) = 

= lim \s{z,)r^'-g'^\s{z2)r^''g''B{tuNk + a-ti) 

11,12-^0 
1 

dxx''~\l - xf^+^^'^^'il -xcos" ^12)^*^ = 
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= lim \siz,)r^'-g'^\s{z2)r^''g''B{t,,Nk + a-t,) 

ii,i2->-0 

B{t2, Nk + a - t2) - ti [ dxx*^~^(l -x)^'=+"~*2"Mn(l -xcos^^ia; 



= lim \s{zi)\-^''g''\s{z2)\~^''g''B{Nk + a,t,)B{Nk + a,t2) 

_ r ^^(1 _ ln(l - x cos^ 6,2) 

Jo X 

Plugging this back in Eq. (93) and taking the limit we get 

Eluj^p{zi)u^p{z2) = 

7 7 ^ riT 

= ^u;^A^iMz2) - ^^ddU^ddU / -(l-x)^^+'^-Mn(l-xcos2M 

Zg -^3 Jo 

and the full expectation value (84) can be expressed using the G-function 

EkU4,p{zi) U^p{z2) = 

1 - - 7 dr 

= uj^^{z,)uj^^{z2) - Tr2dd\.Mz. / da^ — (1 - a;)^'=+'^-Mn(l - xcos^ M = 

^ J-ioo Jo ^ 

= U^^{Z1)U^^{Z2) - 

-^ddUddU r -(1 - xf'^-' ln(l - X cos' eu) ^' (94) 

which is the final answer, valid at finite k. Now we proceed to the analysis of its large 
k behavior. 

5.4- Large k limit 

The integral (94) can be Taylor expanded around cos^ ^12 = 0. Therefore, in the case 
(1) of Section (4.2) of finite Dj{zi,Z2) > the large k behavior of the second term in 
(94) is given by 

^ aai aai dx ^N^~l^ m 2n n ^o[iV°(ir^|l' 2, iVfc) 
- T7:dd\z-^dd\;,,^ / — [1 - xy ln(l - xcos — ~ 



G^^'i{g^^\l,2,...,N,) 



k^ 



which is exactly the behavior predicted by Eq. (71). In order to compute the contact 
terms, we expand (94) for small distances up to the second order in kDi{zi,Z2) 



= ^4>iaih i^l)^<Pia2h (^2) + (95) 

, ^ G^^^f^{g^^\l,2,...,Nk-2,N,-2,N,) ^ 



GS(^^^Hl,2,...,iV, 



^1,^2 7 

ri V i- I I ■ / / \ / , 1 
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see Section 5.6.4 of Ref. [17] for the relevant integrals of the Meijer G-function. Note, 
that the contact terms here have the same structure as in Eq. (76). 

The large k limit value of the ratio of the G-functions in Eq. (95) depends on the 
k dependence of the parameter g = g{k). If — as — oo, then this ratio tends to 
one [20], and we get 

{^4,axaA^l)^<j,a2a2{z2)) g^Q = 

= ^0aiai{Zi)uQa2a2{z2) + (l^Oaiai (^l)l^Oaiai (^l) + I^0a2ai (^l)l^0aia2 (^l) )^^i ,22 ) 

The ratio of the G-functions in (95) is also known if the parameter g scales as 
g ~ A^fc^^, e > 0, see e.g. Ref. [12]. In this case there is an asymptotic expansion of the 
G-function for the large values of the argument z = g'^'' 

GJ:-°(.|6„6„...,6^J = (^^^y''e-^-^^-^.-(l + i^,.-V^^ + (9(.-/^^)), (96) 

where 7 = + jfj^'^n=i^n, Ki ~ A^|. In this case the ratio in Eq. (95) scales as 
^1 g ^ and therefore contact terms vanish at large k 

{uj^{zi)uj^{z2))g^i~^3+. = UJo{Zi)uJoiz2). (97) 

We expect that the same result holds for the standard scaling g Nk [10]. 

6. Conclusions 

In this paper, we considered random Kahler metrics defined through the Bergman metric 
construction by random matrix measures of the eigenvalue type (26). Metric correlators 
for these measures involve nontrivial angular integrals. We were able to compute several 
representations of the two-point functions and in particular perform the analysis of the 
large k asymptotics. We considered in full details the particular example of the Wishart 
measure, where the one and two-point functions can be computed explicitly at finite k. 

The cases we have studied clearly are oversimplified models of random metrics. For 
example, the form of the large k two-point functions (76) at non-coincident points is 
consistent with a 5- function like measure, centered on the background metric uq, when 
viewed as a measure on full space of metrics /C[^(,]. This singular behavior may be 
smoothed to some extent by considering random measures, having ^ log P instead of 
P as the basic (hermitian) random matrix variable [10], which we plan to do in future 
work. However, more interesting and physically relevant models, like the Liouville and 
Mabuchi theories, obviously involve a matrix model measure depending crucially on the 
angular part of the matrix P. In spite of these obvious shortcomings, we believe that our 
present study of the geometric matrix model correlators (5) in the simplest "solvable" 
class of models will be helpful in understanding the more physically relevant cases which 
are presently under investigation. 
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Appendix A. Orthogonal polynomials for the Wishart ensemble 

Here we rederive (94) using the method of orthogonal polynomials. 

It is convenient to scale away the coupling constant g in the Wishart potential, as 
in Eq. (82). Then for the linear potential ^(A) = A the relevant basis of orthogonal 
polynomials (51) is given by the modified Laguerre polynomials [1], satisfying 

Pn(A)p™,(A)e-^A''c/A = /i„(a)5„^, (A.l) 

hn{a) = r(n + l)r(n + 1 + a) 

The integral here converges as long as Rea > — 1, which is true in our case. 
Rescaling \ ^ \/ g in Eq. (46) and using the representation (55) we get 

-—ddU.ddU, lim \s{z,)\~''^\s{z2)\"''' 



' ' ' ' ' (V^lV^2)^-^(^l-^2)Loo^ 



/ A(A)Ai2(A)e~^i^i/^-^^'^2/»e-^-i^«[rfA] ■ Vo\U{Nk 



oo 



—ddl^dd\,, lim |s(^i)|-2ti|s(;z2)r^*' / dxix{'~^ I dx2x'i-^ 



(7/;i7/.2)^^-2(^i-^: 



2 



/7O0 'fc ^ 

■iOO 



d\i \1 / dX2 X2 e-(i+^'i/s)Ai-(i+V2/3)A2 
Jo 

■(PjVfc-2(Al)p7V,,_l(A2) -p7Vfc-l(Ai)pAr^_2(A2)). (A.2) 

Using the Laplace transform formula [1] for the modified Laguerre polynomials 



00 



e-^^AV(A)rfA = ^— TT^r(n + 1 + a) (A.3) 



s 



n+l+a 



we can compute the A-integrals in (A.2) 

) /"OO 
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(PW,-2(Al)piV,.-l(A2) -p7Vfe-l(Al)pjV,~2(A2)) = 



9 



{Nk - l)!(iVfc - 2)! (1 + ^i/g)^^+-{l + ^2/^?)^^+'^ 
Plugging this result back to (A. 2), and rescaling xi^2 9^1,2 we get for the two-point 
function 

EkUJ^p{zi)u^p{z2) = 

-| pioo 

= ^ddU.ddU lim \siz,)r''^g'^\siz2)r'''9'' da g-''-^''^+'^^ 

fOO POO ^fc~ -^ 

■ / dx^x'r' / dX2x'i-'[{l + ^IJi){l + ^lj2)r^'''^''^m n /in(a)-Vol[/(iVfc) = 

7o Jo „=o 



rfa-^9ai,,5ai, lim |.(zor'*^^*^k(^2)r'V^ 

-ioo 1^ h,t2^0 

POO POO 

/ dXix{'"^ C/X2X2'"^(1 +Xi +X2 + XiX2Sin2^12)"^^''-+"^ 

io Jo 



da El uj^p{zi)uj^p{z2), 
which coincides precisely with (92). 
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